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PLASTIC  DEFORMATIONS  OF  A FREE  RING 
UNDER  CONCENTRATED  DYNAMIC  LOADING1 
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A concentrated  time-dependent  force  acts 
on  an  unsupported  thin  ring  along  a diameter#  The 
problem  considered  in  this  paner  is  to  determine  the 
deformations  of  the  ring  when  the  force  magnitudes 
are  such  that  plastic ' strains  large  compared  with 
elastic  strains  occur.  By  neglecting  elastic  strains 
and  assuming  ideally  plastic  behavior,  approximatidns 
to  the  final  deformations  of  the  ring  are  obtained# 

The  analysis  is  developed  for  force  pulses  of  arbi- 
trary shape,  but  numerical  results  are  obtained  only 
in  the  special  case  of  a rectangular  force  pulse.  A 
criterion  is  stated  for  conditions  when  this  type  of 
analysis  can  be  expected  to  provide  satisfactory  re- 
sults# 

l. 

This  analysis  is  an  extension  to  'n  unsupported  ring 
of  ductile  material  of  the  methods  employed  in  These 

references  dealt  with  large  plastic  deformations  of  beams 
subjected  to  dynamic  loading,  under  conditions  such  that  elastic 
deformations  could  be  disregarded# 

The  hypothesis  that  a satisfactory  approximation  to  the 
final  deformation  in  the  ring  subjected  to  a concentrated, 
dynamic  load  (Fig.  2)  can  be  obtained  by  ignoring  elastic 
strains  (and  hence  elastic  vibrations)  leads  to  a "plastic- 
rigid”  type  of  theory#  Briefly,  this  theory  assumes  that  the 
ring  remains  rigid  until  the  moment  developed  in  the  ring# 


1.  The  results  in  this  paper  were  obtained  in  the  course  of 
research  conducted  under  Contract  N7onr-3 5810 'between  the 
- Office  of  Naval  Research  and  Brown  University# 

2#  Research  Associate  in  the  Graduate  Division  of  Applied  Mathe- 
• matics. 

3#  Associate  Professor  of  Engineering# 

*f#  Numbers  in  square  brackets  refer  to  the  bibliography  at  the 
end  of  the  paper. 
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at  some  section  attains  a value  M0*  the  yield  or  limit  moment. 

At  this  section,  unlimited  bending  may  occur  while  the  moment 
equals  Mg.  That  is,  the  tangents  on  either  side  of  this  section 
may  ultimately  make  a finite  angle  with  one  another.  Such  a 
material  is  described  by  Fig.  1.  As  discussed  in  [1],  one  way 
of  computing  MQ  is  by  the  formula 

Mq  = fly  Zp  (l.D 

where  cj  is  the  yield  stress  and  Z is  the  sum  of  the  first 

J r 

moments  of  the  upper  and  lower  halves  of  the  section  area  under 
consideration  about  its  neutral  axis.  When  a section  of  the  ring 
sustains  the  moment  M0,  a plastic  hinge  is  said  to  exist  there, 
and  the  unlimited  bending  that  transpires  corresponds  to  free 
rotation  at  the  hinge  under  constant  motion* 

This  theory  should  yield  useful  approximations  if  the 
total  energy  absorbed  in  plastic  flow  greatly  exceeds  the  elastic 
energy  that  can  be  stored  in  the  ring.  A criterion  expressing 
this  state  of  affairs  is  developed  subsequently. 

It  is  assumed  that  the  middle  line  of  the  ring  has  radius 
r which  is  large  enough  to  permit  the  analysis  to  be  carried  out 
as  if  all  the  material  were  concentrated  at  the  middle  line.  The 
mass  per  unit  length  of  the  middle  line  will  be  denoted  by  m and 
the  time  dependent  load  by  P(t),  The  configuration  to  be  studied 
appears  in  Fig.  2 in  which  small  changes  in  geometry  are  neglected. 

2 . Motion  Immediately  After  _]> i pact 

There  is  uniform  translation  of  the  ring  immediately  after 
the  load  is  applied  and  this  effect  persists  until  the  bending 
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moment  due  to  inertia  forces  attains  the-  value  M0  at  certain 
sections  of  the  ring* 

Let  it  be  assumed  that  these  plastic  regions  first  occur 
at  9 = 0,  7i.  By  symmetry,  the  shear  force  V = 0 at  9 = ft,  and 
V = P/2  at  9 = 0.  The  moment  M(<p)  at  any  cross-section  may  now 
be  obtained  as  indicated  in  Fig.  3„  Here  the  acceleration  q is 
given  by 

P = 2ftrma  , (2.1) 


and 


,9 


Pr 


M(qp)  = NQr(l  - cos  9)  + M0  — — sin  9 + j mar(r  sin  9 - r sin  0)d9 


ft 


2rN0  = 


2 2 
mar  sin  0d0  = 2mar  . 


!0 


These  two  equations  with  (2,1)  give  the  result 


M(cp)  = M [l  - sin  9 + ? sin  9], 


2M, 


2M0 


(2.2) 


If  additional  plastic  regions  occur,  they  will  be  located 
at  stationary  points  of  M(9).  Setting  dH/d9  = 0 gives 


sin  9 = (ft  - 9)cos  9. 


(2.3) 


p p 

From  the  equation  d M/d9^  = Pr/2ft  [(ft  - 9) sin  9+2  cos  9]  it  is 
seen  that  d2M/d92  < 0 only  for  ft/2  « 9 < n,  the  corresponding 
solution  of  (2»3)  being  at  9 = ft#  This  solution  is  already  em- 
ployed. The  solution  of  interest  is 


9 ^ 1.11^  64°  or  sin  9 - (ft  - 9 )cos  9 


(2.4) 


where  d^M/d<^  >0  implying  M(  9)  has  a minimum  there.  Let  the 


1 1 
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load  P be  increased  until  this  minimum  attains  the  value  -MQ  so 


that  (2r2)  gives  the  result 


Pr  _ ^ 6.  38 

% (tt  - Tq) sin  90 


(2.  5) 


which  is  the  collapse  load  parameter  for  the  ring.  When  Pr/MQ 
exceeds  this  value,  the  ring  contains  four  plastic  hinges  at 
<p  = 0,  n,  2it  - £ and  may  be  treated  as  a kinematic  mechanism 
since  the  segments  between  hinges  are  rigid  according  to  the 
basic  hypothesis.  As  in  [l]  it  is  assumed  that  £ = £(t);  in 
special  cases  ^(t)  may  turn  out  to  be  constant, 

3*  Equations  of  Motion  for  the  Deforming  Ring 

Since  the  points  £,  2tc  - f,  correspond  to  stationary  values 
of  M(  cp) , there  is  no  shear  at  these  hinges.  Free  body  diagrams 
of  the  two  segments  joined  at  9 = $ are  given  in  Figs.  4 and  5 
along  with  the  notation  to  be  used  in  the  equations  of  motion. 

The  centers  of  mass  and  moments  of  inertia  with  respect  to  mass 
centers  are  given  by  the  expressions  . 


r = sin 
0 r 2 


ru  = 


n - Z 


I0  = - 2(1  - cos  O j In  ^ |E£[U  - O2  ~ 2(1  + cos  U]j 


Equations  expressing  rate  of  change  of  angular  momentum 
for  each  segment  may  be  written  with  respect  to  transverse  axes 
whose  origins  are  at  the  centers  of  mass  of  the  respective  seg- 
ments. However,  the  origin  for  the  quantities  y0,  y^  is  fixed. 
Although  the  mass  of  each  segment  changes  because  of  the  motion 
of  the  hinge  at  9 = (,  Newton’s  equations  may  be  used  in  the  form 


1 


"summation  of  forces  equals  product  of  mass  and  acceleration  of 
center  of  mass",  and  "summation  of  torques  about  center  of  mass 
equals  product  of  moment  of  inertia  about  center  of  mass  and 
angular  acceleration",  since  the  rate  of  transfer  of  momentum 
from  one  segment  to  the  other  across  the  hinge  due  to  the  chang- 
ing mass  is  zero  and  hence  no  impulsive  forces  are  applied  at 
the  hinge  point. 

The  analysis  is  based  on  the  six  equations  of  motion  for 
the  two  segments  hinged  instantaneously  at  £,  together  with  two 
equations  which  express  continuity  of  velocity  at  the  hinge  point 
5,  It  is  shown  in  Appendix  I that  continuity  of  velocity  ensures 
continuity  of  displacement;  the  acceleration  is  discontinuous, 
in  general,  at  a moving  hinge.  Thus,  eight  equations  are  avail- 

t • 

able  for  the  determination  of  eight  unknowns  y0,  , w0,  u> , £ , 

N0,  N^,  The  first  four  quantities  are  linear  and  angular 

velocities  while  the  last  three  are  axial  forces;  definitions  are 
given  by  Figs,  4 and 

The  equations  of  motion  and  continuity  equations  at  the 
hinge  are  derived  in  Appendix  I making  use  of  Figs,  4 and  5*  The 
six  equations  of  motion  can  be  reduced  to  three  by  eliminating 
N0J  N , Nr.  These  together  with  the  two  continuity  equations 
form  a system  of  five  simultaneous  equations  of  the  form 


ailyo  - ai2yn  + ai3rWo  + 

The  five  equations  (3, 1-3,5) 
values  1,2 ,3, 4,  and  5« 


li4‘ 


'w. 


it 


wo>  wti)  = 0 

3*  2 ) 


are  obtained  by  giving  i in  turn  the 


The  following  notation  wi i 1 be  used: 
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Iffil  (u  - M ) = w 
v -rt  o 7 - 


Pr 


Then  the  coefficients  are  as  follows: 


(3.6) 

(3.7) 


all  = ~ cos  0(5  “ sin  5), 

aX3  = (1  - cos  K)2  - 2 sin  Z(Z  - sin  O, 

al5  = 2 sln  ~ (1  “ cos  Oj 

a22  = (1  + cos  0(tt  - 5 - sin  Oj 


12 

al4 

a21 

a23 


= 0, 


= 0, 


= 0, 


= 0, 


a24  = 

(1  + cos  O2 

- 2 

sin  5(rc  - 

Z - 

sin 

V 

> a25  " " 2 

a31  = 

c» 

a32 

- 

it  - 

a33 

- 

1 - cos  Zf 

a3*t  a 

l + cos 

a35 

— 

- 

2’ 

aifl 

= 

h 

al+2  = 

- 1, 

aM-3 

= 

sin  Z, 

a,.i. 

*T  *T 

= 

- sin  Zf 

ai+5  = 

- ui  cbs  Z , 

a51 

r 

ot 

a52 

rz 

o, 

a53  = 

i - cos  z, 

a5*t 

r 

1 + cos 

0 

a55 

n 

- sin'  Z • 

In  order  to  express  the  accelerations  y0,  y , wQ, 
explicitly  in  terras  of  Z and  p,  let  the  equations  (3.1-5)  he  con- 
sidered as  a system  of  five  algebraic  linear  equations  in  the 

# • • • * * 

four  unknowns  y0,  y^,  rwQ,  rw^  , These  algebraic  equations  will 
be  compatible  and  have  a solution  only  if  the  determinant  of  the 
system  vanishes.  This  condition  is 

I | = 0 (3* ft) 

where  the  elements  in  the  i^*1  row  and  column  are  given  above. 
This  determinant  may  be  split  into  three  determinants  by  breaking 
up  the  last  column.  The  result  is  the  equation 
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Dx(5)  - \ D1(0  - u£  D2(0  = 0 


7 

(3.9) 


where  D^,  D^,  D 2 may  be  reduced  to  third  order  determinants  by 
Ohio's  method  [2]  in  which  form  they  were  used  in  the  numerical 
computations.  These  results  are  listed  in  Appendix  II* 

The  quantity  p + 2 Ku  = D^/D^  } obtainable  from  ( 3 • 9 ) 9 

Is  plotted  against  £ in  Fig.  7. 

The  vanishing  of  (3.9)  means  that  (3.1-5)  has  a unique 
algebraic  solution  in  y*0,  y , rui0,  rufL  which  may  be  obtained  by 
eliminating  the  £ terms  from  (3.*+>5)  and  using  the  resulting  equa- 
tion with  (3.I-3).  This  equation  is 

sin£y0  - sin5yu  + (1  - cos  OruQ  - (1  + cos  Oru^  = 0,  (3. 10) 

Cramer's  rule  will  solve  the  algebraic  system  (3»l~3j  10),  The 
denominator  in  Cramer's  formula  is 


A(0  = !b±3  | 


(3.11) 


"t/tl  • « 

where  the  elements  in  the  i1,  row  and  ;)  * column  are 


bij  " aiJ 


for  i < 3,  1 < b, 


and 


\l  = - \2  ~ sin  b43  = 1 “ °os  £,  b)fi+  = - (1  + cos  £). 


Let  A^(p,0  denote  the  determinant  obtained  from  A(£)  by  replacing 
the  ith  column  of  A(c,)  by  the  negative  of  the  "constant"  terms 
of  (3» l”3j  10))  i. e, , the  terms  not  containing  the  algebraic  un- 
knowns, Then  Cramer's  rule  states  that 


y0  = 


_ A1(p,^) 


aTcT 


••  2 
V =: 


A 0)0 


TC 


A(0 


■>  rio  ' 


A (?) 


A 

ru  = .Jt 

* A(U 


(3.  12) 


I 


1 

j 

. 

I 
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a q.O  - a u,0 
r<“*  " V = a TfT4— ~ 

and  since  A^  and  differ  in  only  one  column,  they  may  be  com- 
bined, giving 

<^(0  = \o±i\  (3.13) 

where  the  elements  in  the  i*'*1  row  and  column  are 

C11  = d “ cos  OU  - sin  O, 

c-^  - O-  - cos  K)(Z  - sin  O, 

.2 


13 

:l4 


(1  - cos  Z)  - 2 sin  Z(Z  - sin  Z ) > 
- 2 sin  £ + ^ ( 1 - cos  S ) , 


C21 

= 

w22 

= (1  + cos 

r Htr 
->  ' **  **■ 

- Z - sin  l) 

? 

C23 

= (1  + cos 

o2  - 

■ 2 sin  j;(m:  - 

S - sin  Z) j 

c24 

= 2 sin 

°31 

= ^ c32 

= 

°33  = 2’ 

34  2’ 

°4l 

= sin  r,  c 

J42 

’ 43 

2 COS  £,  c,^ 

If  the  last  column  of  (3»13)  is  broken  up  so  that  two  determinants 
are  formed,  equation  ( 3 • 13 ) becomes 


wD2U ) = D (Z)  - H d3(U 


(3.14) 


1 ■ 


! 


. 
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where  D^,  may  be  reduced  by  Ohio's  method  to  third  order 
determinants  and  are  listed  in  Appendix  II,  and  DgCO  can  be 
shown  to  be  identical  with  A (O. 

The  derivation  of  the  equations  expressing  continuity 
across  the  hinge  in  Appendix  I shows  that 

(1  - cos  Ow0  + (1  + cos  = 0.  C 3*  X5) 

Consequently,  if  w is  known  u0  and  follow  from  (3*6)  and  (3»15)» 
Then  the  deCormations,  as  shown  in  Appendix  I may  be  obtained  by 
integrating  d©0/dt  = uQ,  d©^/dt  = giving  ©0,  ©^  as  defined 
in  Fig,  6,  The  problem,  therefore,  is  to  solve  the  pair  of  equa- 
tions (3.9)  and  (3.14), 

4#  Solution  for  Rectangular  Force  Pulse 

Numerical  results  for  the  special  case  of  a rectangular 
force  pulse  are  obtained. 

Let  ( = 0 so  that  5 = const.  It  follows  from  (3»9)  that 
p,/2  = D^( O = const,  which  is  necessarily  greater  than 
3,44,  since  all  the  equations  under  consideration  are  valid  only 
when  the  hinge  at  9 = £ is  present,  i,e,  , when  p,  > 6,88  by  (2,5)» 
The  variation  of  $ with  p,  may  be  obtained  from  Fig,  7 on  putting 
K - 0.  Under  the  impact  p-  > 6.83  the  hinge  moves  instantaneously 
to  the  corresponding  position  given  by  Fig,  7*  Let  P = const, 
act  for  a time  T after  which  P = 0,  Let  be  the  corresponding 
hinge  position-  obtained  from  Fig,  7*  Then  w = const,  which, 
from  (3.14)  gives 

D3£0)  - | S «„) 

= “ w 


w 


0 < t < T 


(4.1) 
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where  the  constant  of  integration  is  zero  since  w = 0 for  t = 0, 
Putting  T for  t gives  the  final  value  w of  w under  the  load  P* 


Hj  d3«o>  - \ D (50) 

T = — 


(4.2) 


From  (3.6)  and  (3.15) 


U)  - 


mr- 


1 + cos  ( M 


„0. 


(4.3) 


- 1 - cos  £ M 

It  can  be  seen  from  (4.3)  that  neglecting  co^,  is  equivalent 

p 

to  neglecting  w . Therefore,  by  (4.2),  the  restriction  that  T be 
sufficiently  small  and  decrease  as  P increases  must  be  enforced. 

The  equations  d©Q/dt  = u>0,  d©^/dt  = combined  with  (4.3) 

give 

3 d©.. 

(4.4) 


= -\  (1  + cos 


mt  ;rr  = I'1  - c°s  w** 


(4.5) 


Introducing  (4.1)  into  (4.4,5),  performing  the  integration,  and 
replacing  t by  T gives  the  final  deformations  ©0^>  © ^ under  the 


load. 


3 © 


mL  -OT  _ _ 1 (1  + cos  r ) 

M2  T+  s o t 

o T 


■3  © _ 

mr-5  tcT 


0) 


M T2  = ^ “ cos  50) 

O T 

When  T < t,  P 5 0 and  (3.9,  14)  reduces  to 


(4. 5a) 


(4.  5b) 


Dx(0  = u)CD2(^) 

wD2(0  = D3(^)  or  D^O  = ^ ZD2(f,). 


(4.6) 
(4.  7) 


i 


i 


JU  - «£  -t 


I 
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Dividing  (4,7)  by  (4.6)  gives  the  differential  equation 


11 


D (O 
dw  _ 5 


w 


(4.  8) 


which  may  be  integrated,  (4,2)  serving  as  the  initial  condition, 
giving 


“(5,50)  = exp 

Then,  it  follows  from  (4,6)  that 


i . i oxp(. 

D2(0  ut 


d3(0 

^0 


<n 


(4,9) 


hiiloo 

DX(0 


(‘f.lO) 


Operating  with  the  third  order  determinants  in  Appendix 
II  and  expanding  each  element  in  power  series  in  n/2  - £ it  can 
be  shown  that  = - a(^/2  -5)+...  ,D,>(0=  b + ...  , 

D^( 0 = c + ...  , where  a,  b,  c are  positive  numbers  with 
0 < c/a  <1/2.  Consequently,  f or  £ near  n:/2,  w/T  behaves  like 

p /n 

( Hp/T ) ( it/2  - O and  tends  to  zero  as  K tends  to  u/2.  That  is, 
the  relative  rotation  of  the  two  ring  segments  terminates  when 
£ = x/2  in  view  of  (4.3).  Therefore,  final  deformations  ©0£, 

©^  may  be  obtained  by  integrating  d©Q/dt  s (d©Q/d^)^, 

ft 

d©u/dt  = (d©  /d£)l;  up  to  ( = n/2,  Using  these  latter  relations 
with  (4.4,  5)  gives  the  results 


mr 


3 0 


k/2 


,3  G 


u> 


M. 


0 T 


© 


a£  = aii  _al  _ 1 | (1  t ros  o L (4.  u) 

' — A m ^ J m** 


o t 


>J(p 


TC 


Mo  T2 


^ -f  + i i (1  - cos  ()  L J(, 

o T 1 


n/2 

■i 


w 


u® 


T t, 


(4. 12) 
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where  the  first  terms  on  the  right  of  (4.11,  12)  follow  from 
(4.2)  and  (4.5a,  5b).  Since  behaves  like  (n/2  - £)'1'’‘C//a  when 
Z,  is  near  x/2,  — ^0  as  Z, — ^ -rc/2.  Moreover,  w/t/t^  behaves  like 

(it/2  - ^)2c/a~-1-  so  that  the  improper  integrals  in  (4.11,  12) 
exist. 

Numerical  integration  must  be  used  on  the  above  formulas. 

Since  all  the  integrands  have  infinite  singularities  at  x/2, 

numerical  methods  will  be  used  up  to  Z,  = 1.  5 from  which  point  the 

integrands  are  approximated  to  by  the  terms  in  their  power  series 

up  to  (it/2  - £),  For  this  reason  the  formulas  (4.11,  12)  will  be 

modified.  Put  , . 

10 

(4.13) 


$(£,tp)  = (1  + cos  £)  — 

T t 


u> 


\lf(^,cp)  = (1  **  cos  ^)  . 


(4.14) 


From  the  expansions  given  above  for  D^,  D2 , it  can  be  shown 
that  for  1.  5 <,  Z,  ^x/2 

(1+^-0 [jnUfeSol]  i 1.5  < c,  < it/2 

2 1 (*/2  - I.5)2c/a  (4.15) 


zt  [l  - (f  - 01 


Then  (4.11,  12)  become 
mri  of  _ mr-^  oT 


»<£).!  i L — —1,5  ^ i < tc/2. 

a (n/2  - 1.  5)2c/a  (4. 3.6) 


M 2 
o T 


Mo  f 

p * ^ |_  2 £ i < 

_ l i $>U,<p)d£  - 1 £ ±LLS,J2l\  (_a.  + 2 "_  J_. 

o > 2 a T ■ . ) 


T 


•2c  2c" 


vJqp 


+ 1 


(4. 17) 
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mr 

M, 


•5  6^  3 

£ uf  = mrJ 


S 


itT 


M. 


1 

2 


T 

iM 


\l»(^,9)d^  + -i  & 


mD,*-5  y<£l 


T 


(, 


2c 


JlZ.2  .i«_5 ) 

+ r ’ 


2c /a 


9 


(4.18) 


which  are  the  expressions  for  deformations  used  in  the  computa- 
tions. The  resulting  deformation  parameters  are  plotted  against 
p in  Fig,  8. 


5.  Check  on  Numerical  Results 

When  relative  rotation  of  the  two  ring  segments  ceases, 
the  velocities  y0,  must  be  equal  as  shown  by  (7*12).  Then, 
either  one  may  be  used  to  compute  momentum  of  the  ring.  The  equa- 
tion t*. 

n I 

Pdt  = 2urmy0  = PTtrmy^,  where  t^  is  the  time 

1 


when  the  deformation  terminates,  leads  to  the  formula 


2 v 2 v 

iL  — S3L.  = 2m  UHL. 

^ M0  T 71  T 


(5.D 


since  j1  Pdt  = PT,  From  (3. 12)  and  the  remark  following  (3.14) 
viC 

concerning  Dg( 0 it  follows  that 

2..  D0«)  - 

mr  ,r  _ u 20 

M0  yo  - d2(I)' 


(5.2) 


®iv  -W-jW 

Mo  ^ - 


VO' 


(5.3) 


where  DQ,  DQ  are  derivable  from  A1  (p,0;  D^,  from  A (p,£) 


__  j 

of  (3.12).  These  determinants,  after  reduction  to  3 order,  are 
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listed  in  Appendix  II.  Integration  is  performed  while  \i  is  con- 
stant, the  results  at  t = T being 


^oT  = D°(CP)  " % Po(tP} 

M0  T T>2  CqD 


^ ^ = D*(<||)  - 1 D*w 

Djr?5 


M0  T 


(5.4) 


(5.5) 


For  T < t,  jj.  = 0 and  equations  (5*2,  3)  become 


2 dy  D (0 
mr  o _ o 


TM"  dT  =-^7~T  5 - aU,9} 

im0  at,  t^d0(0 

2 dy  D*(0 

ffi-ar  = -* 5 bu>9>. 

™o  T^D2(C) 


(5.6) 


(5.7) 


Approximations  for  A and  B when  1.  5 <_  £ < m/2  are  obtained  in  the 
same  manner  as  before  with  D0  = c/2  + 2x(x/2  -£)+***  > 

= - c/2  + 2n(x/2  - O + ...  . Then 

A(^cp)  = -JM  -cp )C/a"*1+2it (5  -qp )°/a  V,  1.5  < K < */2 

Ia(S  - 1. 5>o/a t J (5.8) 

= - — - 7 J§(|  -q>)o/a_1-2n(-I  --p)0/a|,  1.5  < 5 < V2. 

Ta(§  - 1. 5)c/at  2 7 J 

2 (5.9) 


Using  these  approximations,  integration  of  (5*6,  7)  up  to  K = n/2 
gives 

ii1’  ^ 


mr  ‘ ~ £o  = mr  ^oT 
V T ‘ T (J(p 

J-.5 


A(^,cp)d^  - joLLl2*slI[1  + ~_(£  -1.5)3(5.10) 

J.  c.  d^C  d 


mr2  y?l  = mr2  YtcT  4 
\ T M0  ' T 


B(c,cp)d^  + i*fLif.i4.^2[i  - |2L-(|  -1.5)1(5.11) 


9 
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Equations  (5*  10,  11)  may  be  introduced  into  (5«l)  as  a check  on 
the  numerical  work  performed#  In  particular,  the  graph  of  the 
right  members  of  (5»D  versus  cp  should  coincide  with  Fig.  7 since 
this  figure  is  essentially  p versus  The  error  shown  there  is 
smallest  for  y^,  as  expected,  because  (5*11)  involves  addition 
only  of  small  numbers.  The  error  in  y0  is  much  larger  since  the 
subtraction  indicated  by  (5*10)  involves  two  large  numbers. 

Error  in  y^  ranges  from  9%  for  large  p to  practically  sero  for 
small  p. 

As  discussed  in  [1]  an  analysis  of  plastic  deformations 
in  which  strains  of  elastic  magnitude  are  neglected  can  be  ex- 
pected to  give  satisfactory  results  when  the  energy  absorbed  in 
plastic  deformation  greatly  exceeds  the  maximum  elastic  strain 
energy  that  could  be  stored  in  the  ring.  This  furnishes  a means 
of  estimating  the  range  of  validity  of  this  type  of  analysis.  The 
maximum  elastic  strain  energy  can  be  estimated  as 

2n  _?  p 2 

P M“  cos  cp  rd{p  M £ 

We  = • 2EI  2 "El* 

The  work  of  plastic  deformation  is  given  approximately  by 


WP  = 2¥0f 

Hence  the  requirement  for  validity  of  the  assumption  that  elastic 
deformations  are  negligible  is 


M?r 

2Mo®0f  »-Hb  • 


(5. 12) 


o 2 

The  dimensionless  deformation  parameter  mr-)9n^/MriT  is  a function 
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of  p-0  only,  as  plotted  in  Fig.  8.  The  inequality  criterion  can 
be  rewritten  as 

4 


rar-^Q 


of  _ 


v2 


^ EXT 


(5. 13) 


4 

,rar  , 

IT  * 


The  period  i of  the  lowest  flexural  mode  of  free  vibration  of  the 
ring  is  given  by 

T = 271  -^2 

Hence  the  inequality  above  can  be  expressed  as 

,3q  « ? 2 

2 


mr^Q.^  n _2 

r = fC,10)  » 20k  2 


Of  - 1 ^ 't  - 0 ll+  T 

~~n  • 


(5. 14) 


This  inequality,  in  conjunction  with  Fig.  8 , provides  a 
convenient  means  of  estimating  the  value  of  the  load  p.Q,  for  a 
given  period  ratio  t/T,  above  which  the  neglect  of  elastic  defor- 
mations in  the  present  analysis  can  be  expected  to  be  Justified. 
Conversely,  for  a given  load  magnitude  1^0  the  above  inequality 
can  be  used  to  estimate  the  upper  limit  of  the  ratio  t/T  above 
which  the  present  results  should  be  valid. 

The  present  results  for  a rectangular  force  pulse  can 
probably  be  used  to  obtain  fair  estimates  of  the  deformations 
which  force  pulses  of  quite  different  shapes  will  produce.  Solu- 
tions of  other  problems  [3],  [4]  have  been  carried  out  for  a 
variety  of  force  pulse  shapes,  e.g,  , a half-sine  pulse,  a trian- 
gular, or  an  exponential  pulse,  as  v/ell  as  for  a rectangular 
pulse.  Comparison  cf  the  results  of  these  solutions  shows  that 
the  curves  of  dimensionless  deformation  (analogous  to  that  plotted 
in  Fig.  8)  for  the  rectangular  pulse  of  duration  I,  and  maximum 


Bll-21 


17 


dimensionless  load  give  reasonably  good  values  of  the  same 

dimensionless  parameter  for  other  force  pulse  shapes,  provided 

in  all  cases  |iQ  stands  for  the  peak  load  and  T stands  for  the 

r® 

time  such  that  n_.T  is  the  total  impulse  \ pdt.  Hence  the  diraen- 

Oo 

sionless  deformations  plotted  in  Fig,  7 can  be  used  tentatively 
to  estimate  the  corresponding  quantity  for  force  pulses  of  general 
shape,  by  interpreting  p,Q  as  the  peak  (dimensionless)  load  and 
T as  the  total  impulse  divided  by  peak  load. 

It  should  also  be  re-emphasized  that  the  present  analysis 
assumes  that  geometry  changes  can  be  neglected  in  the  equations 
of  motion,  so  that  the  results  cannot  be  expected  to  hold  when 
the  distortions  are  very  large.  It  would  be  desirable  to  inves- 
tigate these  and  other  neglected  effects  experimentally  as  well 
as  analytically. 
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Appendix  I 


6.  The  equations  of  motion  for  the  lower  segment  may  be  obtained 
from  Fig.  *f.  These  are 


| - Nj.  sin  q = mr5(y0  + rQ  sin  | u>Q) 

N0  - Nj-  cos  K = - rarOr  - rQ  cos  -^)w0 
2M0  + N0(r  - r0  cos  |)  - N^(r  - r0  cos  |)  - P -2  sin  | 


(6.1) 

(6.2) 


Iowo- 


(6.3) 


After  introducing  into  (6. 1,2,3)  the  values  for  rQ,  I0  from  equa- 
tions (3*0)  it  follows  from  (6.1,2)  that 

NQ  sin  K = £ cos  K - mrS  cos  5y0  + mr2(l  - cos  K - K sin  ()u0 

(6. 4) 

and  from  (6.1)  that 

Njr  sin  S = | - mr^y0  - mr2(l  - cos  (6.?) 

Equation  (3,1)  results  from  combining  (6.3,lf,5). 

Similarly  the  equations  for  the  upper  segment  are  obtained 
from  Fig.  5. 

sin  t,  = mr(u  - O [y ^ + r^  sin  (6.6) 

N,r  cos  K - Nu  = mr(ai  - 0(r  - r^  cos  ■-  (6.7) 

-2M0  + N^(r  - r^  cos  --“—£)  + N^(r  - r^  cos  JL_n_I)=lnw^.  (6.8) 
Combining  (6.6,  7),  after  introducing  r,(,  1^,  from  (3.0)  gives 


N^sin  5 = mr(  71  - ^)cos^y7i;  + mr2[l  + cos  £ -(it  - Osin  K (6.  9) 

and  substituting  this  equation  and  (6,6)  into  (6.8)  gives  equation 
(3.2). 
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7.  The  Equations  for  displacement,  velocity  and  acceleration  at 
the  moving  hinge  may  be  obtained  in  the  following  manner,, 

Let  © = 9(9, t)  denote  the  small  angle  through  which  each 
element  rd9  of  the  ring  has  rotated*  Deformation  due  to  tensile 
stresses  will  be  neglected  so  that  this  rotation  represents  the 
deformation  of  each  ring  element.  It  follows  from  Fig.  9 that 

dy  = rd9  sin(©  + 9)^  r©  cos  9d9  + r sin  9d9  (7»l) 

dx  = rd9  cos(9  + 9) ^ r cos  939  - r©  sin  9d9.  (7*2) 

The  existence  of  the  plastic  hinge  at  9 = ^ implies  dif- 
ferent angular  velocities  for  the  two  ring  segments  joined  at  the 
hinge.  In  general,  then,  the  deformation  function  ©(9,t)  will 
be  different  for  the  two  ring  segments.  Let  ©0  = ©(9,t),  9 < £ 
and  ©^  = ©(9 , t) , ( < 9 so  that  9©0/3t  = u>Q(t),  90^/at  = w^'t). 

Note  that  ©0  and  ©^  may  depend  on  9,  but  w0,  which  describe 
the  rigid  motion  of  the  lower  and  upper  segments  respectively  of 
the  ring,  must  depend  on  t alone.  That  these  ring  segments  are 
rigid  follows  from  the  rigid-plastic  hypothesis.  For,  the  moment 
M(9 ) in  the  portion  of  the  ring  over  which  the  hinge  at  £ has 
passed  is  less  than  the  limit  moment  M0,  The  alternative  is  that 
M(9)  = M0  and  then  V = dM/rd9  = 0.  But  V ^ 0 because  of  the  pres- 
ence of  d'Alembert  fobces  proporational  to  the  local  acceleration. 

The  deformation  (relative  rotation  of  the  ring  segments) 
generated  at  the  moving  hinge  appears  in  the  form  of  an  increase 
in  curvature  since  the  hinge  remains  at  no  point  for  a finite 
time.  This  increase  in  curvature  is  given  by 


9 

nr  dt  - -3T  dt  **  " “b 


r£dt 


(7.3) 
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Moreover,  as  long  a3  the  hinge  moves  no  "kinks"  appear  so  that 


e0(«f,t)  = ,t). 


(7.4) 


It  follows  from  (7.1)  that 


$ 


„9 


y = yQ  + r | ©0  cos  cpdcp  + r | sin  cpd9,  9 < £.  (7.5) 

Jo  Jo 

Differentiation  with  respect  to  time  is  denoted  by  dots.  From 

(7.5)  <? 


y = yQ  + rwQ  : cos  cpdcp,  9 < *, 

Jo 

9 


♦ • •• 


(7.6) 


(7.7) 


y = y0  + rojQ  j cos  9d9,  9 < 

JO 

Remark:  The  approximation  made  in  (7.1,  2)  is  equivalent  to 

neglecting  u)2,  a restriction  mentioned  previously  under  (4.3), 
for  if  (7.7)  were  obtained  without  making  this  approximation,  an 
additional  term  containing  the  factor  would  appear. 

When  K < 9,  it  again  follows  from  (7.1)  that 


l 

1 r 


9 

(V 


y = y0  + r 


©0  cos  9d9  + r 


00 


©^  cos  od9  + r 


sin  9d9,  ^<9  (7.  8) 


0 


K 

I' 


y = y0  + r©0U",tK  cos  C - r©7l(^+, t)^  cos  K + rwQ  | cos  9d9 

JO 

? 


+ ru^  ■ cos  9d9,  £ < 9, 


(7.9) 


This  last  equation  simplifies  on  using  (7.4)  or,  if  the  hinge 


1 


1 

I 


I 
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velocitv  is  zero,  deleting  the  terms  containing  K so  that 

rf  f 

y = y0  + rwQ  j cos  <pd9  + ru^ 


£ 


»•  *• 


Jo 


k 


y = y0  + ru  i cos  qpdcp  + ri^ 


cos  9dq>,  £ < 9 (7*10) 

cos  9d9  + r cos  C(wQ-  *%)£,  £ < 9* 

(7.11) 


Now  let  9—*  C in  (7.?,  6.  7)  and  9 — > in  (7.«,  10,  11).  The 

results  are  that  y”  = y+,  y“  = y+,  but  y+  - y~  = r cos 
That  is,  vertical  components  of  displacements  and  velocities  are 
equal  at  the  hinge,  but  accelerations  are  discontinuous,  and  they 
jump  by  the  amount  r cos  S(u0  - «*>K)^. 

Putting  9 = 7i,  y = y^  in  (7.10,  11)  yields 


M •• 

y0  - 7* 


y*  = y0 

+ rwQ  sin  K 


+ ru0  sin  5 
- rw^  sin  5 


ru  sin  £ (7.12) 

Tl 

9 

r cos  ^(w^  - = 0.  (7.13) 


A similar  analysis  applies  to  the  horizontal  components. 
Prom  (7.2) 


x = r 


i?  r? 

cos  9d9  - r I ©0  sin  9d9,  9 < S 


(7.  lb) 


x = r 


9 

V 

n r 

cos  9d9  - r 

©n  sin  9d9  - r 

Jo  J 

0 Jl 

©n  sin  9d9,  t,  < 9.  (7.1?) 


Differentiating  these  last  two  expressions  yields  the  results 
that  horizontal  components  of  displacement  and  velocity  are  con- 
tinuous across  the  hinge  but  acceleration  jumps  by  the  amount 
r sin  5(  - uQ) 

Putting  9 = ■*,  x = x^*0  into  the  differentiated  forms 
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of  (7.14,  1?)  gives 

(1  - cos  Oru0  + (1  + cos  Oru^  = 0 (7«  16) 

• • • 

(1  - cos  Orw0  + (1  + COS  Oruj  - sin  ^rCw^  - w0K  = 0,  (7.17) 

Remark:  (i)  Equations  (7*12,  16)  may  be  obtained  by  requiring 

velocities  to  be  equal  at  the  hinge  and  (7.13,  17)  may  then  be 
obtained  by  differentiation. 

(ii)  The  quantities  0Q,  defined  by  = 3©0/5t, 
wit  = d®n/9t  are  required.  These  deformations  will  be  obtained 
only  at  cp  = 0,  as  shown  in  Fig,  6.  Since  the  hinge  at  5 never 
reaches  9 = 0,  it,  the  desired  deformation  may  be  calculated  from 
d©0/dt  = 4-,,  dO^/dt  = because  of  the  rigid  property  of  the 
ring  in  the  vicinity  of  these  points. 
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Appendix  II 


8.  Below  are  forms  of  the  various  determinants  as  used  in  the 
computations  where  the  elements  in  the  i^  row  and  column  are 
denoted  by  a^. 


where 

du 


D0(O  = 2 sin£(l  + cos  O I di;j  I (8.1) 


l21 


13 


l22 


23 


31 


1,  dl2  = sin  r,(i-cos  O - 2(l+cos  £)(s-sin  O 
(1  - cos  O2  - 2 sin  - sin  O, 

(1  + cos  £)(n  - K - sin  O 

- ^_|1^^[2(  1-cos  0(^-0  sin  O - sin  £(l+cos  Ol 
0,  d^^  = tt  - £ + sin  S,  d^  = 2(1  - cos  O 


D0(O  = (1  + cos  O lei;j  1 

(8.2) 

where 

eU  = 1 - cos  e2l  = °»  e3l  = eij  = dij 

for  j = 2,3 

D^(  O = 2(1  + cos  Osin  Clf^l 

(8.3) 

where 


fll  = " f21  = 1?  f12  = - (1  - °os  ~ sin  O,  f“13  = d13> 

f 22  = “ sin  Ol  + cos  Q + 2(1  - cos  0(TC“S-sin  0,  f23  = d23, 
f^]_  = °>  f32  = - + sin  0»  -33  = 2(1  - cos  Q 

Dr(0  = (1  + cos  O i gij  i 


(8.4) 
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where 


2? 


gH  = ! - cos  5,  gp-i  = 0,  g,T  = 1,  for  j = 2,3 


>21 


31 


’±1  “ ij 


where 
h 


D^U)  = 2(1  + cos  £ ) | | (8.5) 

^ - (1  - cos  OU  - sin  Oj  h12  = sin  £(1  - cos  O-'+CS-sinS  ) 


*13 


- h2^  = 1,  h21  = (1  + cos  0(*  - K - sin  O 


h22  = 2(1  - cos  0(tt  -5”sin  O - sin  S(1  + cos  S), 


h31 

where 

kiJ 

where 


= it(l  - cos  ?),  h^2  = - 2£(1  - cos  O,  h-^  = 0 

Dx(0  = sin  Slkjj  i 

= h, . for  J = 1,2  and  k = k =1,  k_^  = 0 

ij  13  33  23 

D2(  O = sin  1 


for  j = 1,2  and  *1— i — JlifLil j; 


(8.6) 


(8.7) 


"13 

= hid 

for 

^23 

= o, 

^33 

where 

mll 

= hll» 

m12 

m21 

= h21» 

m22 

m3l 

= V 

m32 

where 

nij 

= ffiij 

for 

1 - cos  £ 


D 


3(0  =2  sin  S|m1;.  ( 


13 


(3.8) 

: " m23: 


D ( 0 = ( 1 - cos  0 |n^  | 


= m.j  ^ for  ;}  = 1,2  ana  = 1,  m 


(8.9) 
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Asst.  Chief  for  Military 
Construction 
Department  of  the  Army 
Bldg.  T-3,  Gravelly  Point 
Washington  25,  D.  C. 

Attn:  Structures  Branch 
(M.  F.  Carey) 

Protective  Construction 
Branch  (I.  0.  Thornley) 

Office  of  the  Chief  of  Engineers 
Asst.  Chief  for  Military 
Operations 

Department  of  the  Army 
Bldg.  T-7,  Gravelly  Point 
Washington  25,  D.  C. 

Attn:  Structures  Development 
3ranch  (W.F.  Woollard) 

U.S.  Army  Waterways  Experiment 
Station 
P.  0.  Box  631 
Halls  Ferry  Road 
Vicksburg,  Mississippi 
Attn:  Col.  H.  J.  Skidmore 

The  Commanding  General 
Sandia  Base,  P.  0.  Box  5100 
Albuquerque,  New  Mexico 
Attn:  Col.  Canterbury 

Operations  Research  Officer 
Department  of  the  Army 
Ft.  Lesley  J.  McNair 
Washington  25,  D.  C. 

Attn.  Howard  Brackney 

Office  of  Chief  of  Ordnance 
Office  of  Ordnance  Research 
Department  of  the  Army 
The  Fontagon  Annex  7!'2 
Washington  25,  D.  C. 

Attn:  ORDTB-PS 

Ballistics  Research  Laboratory 
Aberdeen  Froving  Ground 
Aberdeen,  Maryland 
Attn:  Dr.  C.  V/.  Lampson 
c)  Navy 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  OP-31 
OP -3^3 


(1) 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
’Washington  25,  D.  C. 

Attn:  Director  of  Research 
Code  423 
Code  442 
Code  421 


(1) 


(1) 


(1) 


(2) 

(1) 

(1) 

(1) 


(1) 


(1) 


Director,  David  Taylor  Model  Basin 
Department  of  the  Navy 
Washington  7,  D.  C. 

Attn:  Code  720,  Structures 

Division  (1) 

Code  740,  Mi-Speed 

Dynamics  Div.  (1) 

Commanding  Officer 

Underwater  Explosions  Research  Div. 

Code  290 

Norfolk  Naval  Shipyard 
Portsmouth,  Virginia  (1) 


Commander 

Portsmouth  Naval  Shipyard 
Portsmouth,  N.  H. 

Attn:  Design  Division 

Director,  Materials  Laboratory 
New  York  Naval  Shipyard 
Brooklyn  1,  New  York 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D,  C. 

Attn:  Ad -3,  Technical  Library 
Rec,  P.  H,  Girouard 

Naval  Ordnance  Laboratory 
White  Oak,  Maryland 
RFD  1,  Silver  Spring,  Maryland 
Attn:  Mechanics  Division 
Explosive  Division 
Mech.  Evaluation  Div. 

Commander 
U.S.  Naval  Ordnance  Test  Station 
Inyokern,  California 
Post  Office  - China  Lake,  Calif. 
(1)  Attn:  Scientific  Officer 


(1) 


(1) 


(1) 


(1) 

(1) 


(1) 

(1) 

(1) 


(1) 


(1) 

(1) 


Naval  Ordnance  Test  Station 
Underwater  Ordnancb  Division 
Pasadena,  California 
Attn:  Structures  Division  (1) 


! ; 


1 
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Navy  (cont.) 

Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  TD-4l,  Technical  Library 

(1) 

Chief.  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  P-314  (1) 

Code  C-313  (1) 

Officer  in  Charge 
Naval  Civil  Engr.  Research  <?* 
Evaluation  Laboratory 
Naval  Station 

Port  Hueneme,  California  *.1) 
Superintendent 

U.S.  Naval  Post  Graduate  School 
Annapolis,  Maryland  (1) 


d ) Ai_r  Pj3rce_s 

Commanding  General 
U.S.  Air  Force 
The  Pentagon 
Washington  25,  D.  C. 

Attn;  Res.<?-  Development  Div.(l) 

Deputy  Chief  of  Staff,  Operations 
Air  Targets  Division 
Headquarters,  U.S.  Air  Force 
Washing ton  25,  D.  C. 

Attn:  AFOIN-T/PV  (1) 

Office  of  Air  Research 
Wright -Patterson  Air  Force  Base 
Dayton,  Ohio 

Attn:  Chief,  Applied  Mechanics 

Group  (1) 

e ) Other.  Government  Agenc ies 

U.S.  Atomic  F.nergy  Commission 
Division  of  Research 
Washington,  D.  C.  (1) 

Director,  National  Bureau  of 
Standards 
Washington  25,  D.  C. 

Attn:  Dr.  W.II.  Romberg  (1) 


S_ujjplejeen_tar_y_  D i s t. r_ ibu; ti on. JL i st 
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Unclassified  Classified 

P.epjqrts Reports 

Professor  Lynn  Beedle 
Fritz  Engineering  Laboratory 
Lehigh  University 

Bethlehem,  Pennsylvania  1 

Professor  R.L.  Bisplinghoff 
Dept,  of  Aeronautical  Engineering 
Massachusetts  Institute  of  ,|''echnology 

Cambridge  39,  Massachusetts  1 1 

Professor  Hans  Bleich 
Dept,  of  Civil  Engineering 
Columbia  University 
Broadway  at  117th  St. 

New  York  27,  New  York 


1 
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Unclassified 
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Classi 
Uejsor  t 


Professor  B. A.  Boley 

Dept,  of  Aeronautical  Engineering 

Ohio  State  University 

Columbus,  Ohio  1 

Professor  G.F,  Carrier 
309  Fierce  Hall 
Harvard  University 

Cambridge,  Massachusetts  1 

Professor  R.J.  Dolan 

Dept,  of  Theoretical  f-  Applied 

Mechanics 

University  of  Illinois 

Urbana,  Illinois  1 

Professor  Lloyd  Donnell 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Technology  Center 

Chicago  16,  Illinois  1 

Professor  A.C.  Er ingen 
Illinois  Institute  of  Technology 
Department  of  Mechanics 
Technology  Center 


Chicago  16,  Illinois  1 

Professor  B.  Fried 

Dept,  of  Mechanical  Engineering 

Washington  State  College 

Pullman,  Washington  1 

Mr.  Martin  Coland 

Midv/est  Research  Institute 

4049  Pennsylvania  Avenue 

Tansas  City  B,  Missouri  1 

Dr.  J.N.  Goodier 
School  of  Engineering 
Stanford  University 

Stanford,  California  1 

Professor  1.!!.  TTe rmes 

College  of  Engineering 

University  of  Santa  Clara 

Santa  Clara,  California  1 

Professor  R.J.  Hansen 
Dept,  of  Civil  Sanitary  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge  3%  Massachusetts  1 


1 


1 ' 


1 
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Addressee  Unclassified 

Reports  

Professor  M.  Hetcnyi 
Walter  P.  Murphy  Lrofessor 
Northwestern  University 
Evanston,  Illinois  1 

Dr.  N.J.  Hoff,  Head 
Department  of  Aeronautical 
Engineering  <?:  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
Brooklyn  2,  New  York  1 

Dr.  J.H.  Hollomon 

General  Electric  Research  Laboratories 
1 River  Road 

Schenectady,  New  York.  1 

Dr.  W.H.  Koppmann 

Department  of  Applied  Mechanics 

Johns  Hopkins  University 

Baltimore,  Maryland  1 

Professor  L.S.  Jacobsen 
Department  of  Mechanical  Engineering 
Stanford  University 

Stanford,  California  1 

Professor  J.  Kempner 

Department  of  Aeronautical  Engineering 
and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
99  Livingston  Street 

Brooklyn  2,  New  York  1 

Professor  George  Lee 

Department  of  Aeronautical  Engineering 
Renssalaer  Polytechnic  Institute 
Troy,  New  York  1 

Professor  Paul  Lieber 
Department  of  Aeronautical  Engineering 
Renssalaer  Polytechnic  Institute 
Troy,  New  York  1 

Professor  Glen  Murphy,  Head 
Department  of  Theoretical  <?: 

Applied  Mechanics 
Iowa  State  College 

Ames,  Iowa  1 

Professor  N.M.  Nev/mark 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  1 
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Unclassified 
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Professor  Jesse  Ormondroyd 

University  of  Michigan 

Ann  Arbor,  Michigan  1 

Dr.  V/.  Osgood 

Armour  Research  Institute 

Technology  Center 

Chicago,  Illinois  1 

Dr.  R.P.  Petersen,  Director 
Applied  Physics  Division 
Sandia  Laboratory 

Albuquerque,  Nev;  Mexico  1 

Dr.  A.  Phillips 
School  of  Engineering 
Stanford  University 

Stanford,  California  1 


Dr.  W.  Prager 

Graduate  Division  of  Applied  Mathematics 
Brown  University 

Providence  12,  R.  I.  1 

Dr.  S.  Raynor 

Armour  Research  Foundation 
Illinois  Institute  of  Technology 
Chicago,  Illinois  1 

Professor  E.  Eeissner 
Department  of  Mathematics 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  1 

Professor  M.A.  Sadov/sky 
Illinois  Institute  of  Technology 
Technology  Center 

Chicago  16,  Illinois  1 

Professor  V.L.  Salerno 
Department  of  Aeronautical  Engineering 
Renssalaer  Polytechnic  Institute 
Troy,  New  York.  1 

Professor  M.G.  Salvador! 

Department  of  Civil  Engineering 

Columbia  University 

Broadway  at  117th  Street 

New  York  2T,  New  York  1 

Professor  J.E.  Stallmeyer 

Talbot  Laboratory 

Department  of  Civil  Engineering 

University  of  Illinois 

Urbana,  Illinois  1 
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Unclassified  Classified 

Reports Reports 


Professor  E.  Sternberg 
Illinois  Institute  of  Technology- 
Technology.  Center 

Chicago  16,  Illinois  1 

Professor  R.  G,  Sturm 
Purdue  University 

Lafayette,  Indiana  1 

Professor  F.  K.  Teichmann 

Department  of  Aeronautical  Engineering 

New  York  University 

University  Heights,  Bronx 

New  York,  N*  Y.  1 

Professor  C.  T.  V?ahg 

Department  of  Aeronautical  Engineering 

New  Yoik  University 

University  Heights,  Bronx 

New  York,  N.  Y.  1 
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